Anomalous Scaling and Refined Similarity of an Active Scalar in a Model of 

Homogeneous Turbulent Convection 
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Anomalous scaling in the statistics of an active scalar in homogeneous turbulent convection is 
studied using a dynamical shell model. We extend refined similarity ideas for homogeneous and 
isotropic turbulence to homogeneous turbulent convection and attribute the origin of the anomalous 
scaling to variations of the entropy transfer rate. We verify the consequences and thus the validity 
of our hypothesis by showing that the conditional statistics of the active scalar and the velocity 
at fixed values of entropy transfer rate are not anomalous but have simple scaling with exponents 
given by dimensional considerations, and that the intermittency corrections are given by the scaling 
exponents of the moments of the entropy transfer rate. 
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Since the work of Kolmogorov in 1941 (K41) 0, much 
effort has been devoted to the study of the possible uni- 
versal statistics of fluid turbulence in the inertial range, 
the range of length scales that are smaller than those of 
energy input and larger than those directly affected by 
molecular dissipation. A major challenge is to under- 
stand, from first principles, the origin of anomalous scal- 
ing, which is the deviation of the velocity scaling behav- 
ior from those predicted by dimensional considerations in 
K41. One important idea proposed by Kolmogorov in his 
refined theory 2], which we refer to as Kolmogorov's re- 
fined similarity idea, attributes the origin of anomalous 
scaling of the velocity to the variations of local energy 
dissipation rate. Kraichnan Q later pointed out that 
the local energy dissipation rate is not an inertial-range 
quantity and proposed to attribute the origin of anoma- 
lous scaling of the velocity instead to the variations of 
the local energy transfer rate and we shall refer this as 
Kraichnan's refined similarity idea. 

Similar problems of anomalous scaling can be posed 
for a scalar field advected by a turbulent velocity field. 
A passive scalar leaves the velocity statistics intact while 
an active scalar couples with the velocity and influences 
its statistics. The problem of anomalous scaling of pas- 
sive scalars is linear and has been recently understood 
in terms of statistically preserved structures 0. On the 
other hand, the nonlinear problem of anomalous scaling 
of active scalars, like that of velocity, remains unsolved. 
A common example of an active scalar is temperature 
in turbulent convection in which temperature variations 
drive the flow. Turbulent convection is often investi- 
gated experimentally in Rayleigh-Benard convection cells 
heated from below and cooled on top (see, e.g., 0,0,0 f° r 
a review). Such confined convective flows are highly in- 
homogeneous with thermal and viscous boundary layers 
near the top and the bottom of the cell. Moreover, coher- 
ent structures, known as plumes, could affect the scaling 
properties For the purpose of studying anomalous 



scaling of an active scalar, it would thus be more desir- 
able to study homogeneous turbulent convection and in 
the absence of coherent structures. 

Homogeneous turbulent convection has been pro- 
posed |9j as a convective flow in a box, with periodic 
boundary conditions, driven by a constant temperature 
gradient along the vertical direction. In Boussinesq ap- 
proximation [lo| , the equations of motion read [hO ] : 



du 

at 

09 

at 



u ■ Vu 



- Vp + u + agOz 



u-V9 = kV 2 9 + (3u z 



(1) 
(2) 



with V • u = 0. Here, u is the velocity, p is the pressure 
divided by the density, 9 = T — (To — f3z) is the deviation 
of temperature T from a linear profile of constant tem- 
perature gradient of —f3, Tq, a, v, and k are respectively 
the mean temperature, the volume expansion coefficient, 
kinematic viscosity and thermal diffusivity of the fluid, g 
is the acceleration due to gravity, and z is a unit vector in 
the vertical direction. The Bolgiano length [12j , given by 
Lb = e 5 / 4 x~ 3 / 4 (ai7)™ 3 / 2 , where e and x are respectively 
the average energy and thermal dissipation rates, is an 
estimate of a length scale above which buoyant forces are 
dominant. Numerical studies 0, [3] revealed that Lb is 
of the order of the size of the periodic box. As a result, 
temperature is not active in the intermediate scales. In- 
deed the small-scale isotropic fluctuations were found [l3| 
to have scaling close to that of K41. 

On the other hand, a dynamical shell model for homo- 
geneous turbulent convection driven by a temperature 
gradient has also been proposed [l4| . Shell model is con- 
structed in a discretized Fourier space with k n — k§h n , 
n = 0,1, . . . , N — 1, being the wavenumber in the nth 
shell, and h and kg are customarily taken to be 2 and 1 
respectively. Shell models for homogeneous and isotropic 
turbulence have been studied extensively and proved to 
be successful in reproducing the scaling properties ob- 



served in experiments 15[. In this shell model for ho- 
mogeneous turbulent convection, real variables u n and 
9 n representing the Fourier transforms of u and 9 satisfy 
the evolution equations: 

^- + vk 2 n u n = ak n (u 2 l _ 1 - hu n u n+ i) 
at 

+ bk n { - hu 2 n+l ) + ug9 n (3) 

d9 

—^- + nk 2 l 9 n = ak„(u n -i9 n ^i - hu n n+ i) 
at 

+ bk n {u n 9 n -i - hu n+ i9 n+1 ) + 0u n (4) 



where a, 6, a, and b are positive parameters. It was 
reported that the scaling behavior depends on the ratio 
b/a 0: close to Bolgiano-Obukhov (BO) scaling [HEl] 



(u n ~ k n 3 ^ 5 , 9 n ~ k n 1/0 ) for b/a large and close to K41 



-1/5 



9 n ~ k n 1 ^ 3 ) for b/a smaller than 



scaling Q (u n ~ fc„ 1/3 , 
about 0.4. In this Letter, we show that buoyant forces 
are important when b/a is large and insignificant when 
b/a is small in this shell model of homogeneous turbulent 
convection and focus on b/a large to study anomalous 
scaling of an active scalar. We extend refined similarity 
ideas for homogeneous and isotropic turbulence to ho- 
mogeneous turbulent convection. Specifically, we extend 
Kraichnan's refined similarity idea and attribute the ori- 
gin of the anomalous scaling in homogeneous turbulent 
convection to variations of the entropy transfer rate. Us- 
ing numerical simulations of the model, we verify the 
consequences and thus the validity of our hypothesis. 
Multiply Eq. §5§ by u n , we get the energy budget: 



dE„ 



F u {k n ) - F u (k n+ i) - vk n u n + agu 



(5) 



dt 

where E n = u^/2 is the energy in the nth shell, F u (k n ) = 
k n (au n -i+bu n )u n ^iu n is the rate of energy transfer from 
(n— l)th to nth shell, vk^u^ is the rate of energy dissipa- 
tion in the nth shell due to viscosity, and agu n 9 n is the 
power injected into the nth shell by the buoyant forces. 
It is thus reasonable to take buoyancy to be significant 
in the nth shell if 



ag\(u n 9 n ) \ > e 



(G) 



where (. . .) is an average over time. Note that for both 
K41 and BO scaling, ag(u n 9 n ) = e at k n = 1 /Lb- We 
find that Eq. ([6]) is satisfied for most of the shells and 
thus 9 n is active only for b/a is greater than about 2 (see 
Fig. H]), and that this change coincides with the reported 
change in scaling behavior discussed above. 

As we are interested in the case of an active scalar, we 
focus on b/a large and study the velocity and tempera- 
ture structure functions, S p (k n ) and R p (k n ): 

"<* • R p (k n ) = (\9 n \P) ~ 



s p {k n ) = (Kn 



(7) 



The scaling exponents (, p and £ p do not depend on the 
values of the various parameters as long as b/a is larger 
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FIG. 1: The logarithm of ag\{u n 0„)\/e for different shells for 
a = 0.01, & = 1, /3 = 1, v = 5 X 10" 17 , s = 5x 10" 15 and 
N = 32 (circles), and a = 10,6 = 1, f3 = 100, v = k = 10~ 8 , 
and JV = 30 (squares) (the datapoints for n — 28 and 29 
are not shown here as they are too small. For both cases, 
a = b — 1 and ag = 1. 



than 2. The results reported below are obtained using 
a = 0.01, 6=1, /3 = 1, 5 = 6 = 1, ag = l, i/ = 5x 10~ 17 , 
k = 5 x 10~ 15 and N — 32. As can be seen in Fig. [2l 
both £ p and £ p deviate respectively from the BO values of 
3p/5 and p/5 obtained from dimensional considerations. 
Thus the active temperature and the velocity in homoge- 
neous turbulent convection have anomalous scaling. We 
study also the case where the temperature is driven by a 
large-scale random forcing instead of an imposed linear 
gradient. In this case, the j3u n term in Eq. ^ is re- 
placed by a random noise acting only in shell n = 0. We 
find exactly the same scaling exponents, supporting the 
universality of scaling of an active scalar upon different 
forcing mechanisms 17, 18j |. 
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FIG. 2: Deviation of the scaling exponents from the BO val- 
ues: — 3p/5 (circles) and£ p — p/5 (squares). The solid lines 
are the results of Eqs. (jT5J and ([To]) . 



It was suggested [19( that when buoyancy is dominant, 
the scaling behavior of velocity and temperature spectra 
is governed by an entropy cascade of constant entropy 
flux. In Bousinessq approximation, the entropy is pro- 
portional to the volume integral of temperature fluctua- 
tions. Entropy in the nth shell is, therefore, defined as 
S n = 9^/2. By studying the entropy budget obtained 
from Eq. |4|) upon mulitplication by 9 n : 

-0± = F e (k n ) - F e (k n+1 ) - nklel + [3u n 9 n (8) 

we get the rate of entropy transfer or entropy flux from 
(n — l)th to nth shell as: 
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F g (k n ) = k n (au n _ 1 + bu r , 



(9) 



In the stationary state and for intermediate scales where 
scaling is observed, both /3{u n n ) and nk 2 i (8 2 l } are neg- 
ligible such that there is indeed a constant entropy flux 

with <F e (AvO}=X = 

We propose that when buoyancy is significant, 

u n = M^9) 2/5 \Fe(k n )\ 1/5 k0 5 (10) 

9 n = <t>o{ag)- 1 '*\F e (K)? l *K y * (11) 

where <f> u and <j)g are dimensionless random variables that 
are independent of k n and statistically independent of 
Fg(k n ). The absolute signs are taken because the entropy 
flux Fg(k n ), unlike can assume both positive and neg- 
ative values. Equations. (fit)]) and (111|) are an extension of 
Kraichnan's refined similarity idea to homogeneous tur- 
bulent convection. With Eqs. (|10|) and (TTTj) , we attribute 
the anomalous scaling behavior of the active tempera- 
ture and the velocity to the shell-to-shell variations of 
the entropy transfer rate. An immediate consequence is 
that the conditional velocity and temperature structure 
functions at a certain prescribed value x of the entropy 
transfer rate are given by 



{\u n \ p \F e (k n ) = x) = 
(K\ p \F e (k n ) = x) = 



kn Cp (12) 



<</>£> (asr f z* fc n f ~ k n Cp (13) 



and hence would have simple scaling with BO exponents 
of Cp = 3p/5 and £* = p/5 respectively. We evaluate the 
conditional velocity and temperature structure functions 
at different values of x and confirm that Q and £* are 
independent of x, and in good agreement with 3p/5 and 
p/5 respectively as shown in Fig. [3] 



Let (\Fg{k n )\ p ) ~ k n Tp , then Eqs. (JTUJ) and (JTTJ) imply 

C P = 3p/5 + r p/5 ; £ p = p/5 + T 2p/5 (14) 

showing that the intermittency corrections, which are the 
deviations of the scaling exponents from the BO values, 



FIG. 3: The scaling exponents Q (diamonds) and Q (tri- 
angles) of the conditional velocity and temperature structure 
functions. They are in good agreement with the BO values of 
3p/5 and p/5 (dashed lines). 



are given by the scaling exponents of the moments of the 
entropy transfer rate. As the power of Fg in R p is twice 
that in S p , this explains why the anomaly is larger for £ p 
than for ( p (see Fig. [5]). We evaluate r p numerically and 
check Eq. (fT4"l) in Fig.0] Good agreement is again found. 




FIG. 4: Comparison of £ p (circles) and £ p (squares) with 
the theoretical predictions of 3p/5 + t p / 5 (stars) and p/5 + 
T 2 P /5 (crosses) using the numerical results of t p . 

Next, we show that the intermittency corrections, as 
given by t p , can be obtained by suitably modifying 
the results of the scaling exponents of the moments 
of the local thermal dissipation rate found in experi- 
ments [20[ ■ In Ref . [2(j , the statistics of the local ther- 
mal dissipation rate, estimated by % T , have been stud- 
ied in the central region of turbulent Rayleigh-Benard 
convection. It was found that the moments of \ T = 
((Wc)t)" 1 Jt +T n{dT/dt') 2 dt', where {u 2 c ) is the mean 
square velocity fluctuations at the centre, satisfy a hi- 



4 



erarchical structure of the She-Leveque form [21], and 
that their scaling exponents fi p , defined by (x%) ~ T>lp > 
can be well described by [i v = c(l — (3V) — Xp with c = 1, 
(3 X = 2/3 and A = 1/3. The parameter A is the scaling ex- 
ponent of lim p _ >00 (x^ +1 )/(x^), which was estimated [2(| 
as the ratio of the maximum thermal dissipation divided 
by a time t r at the scale r = y/ (u|)t. Taking t r as 
r/u r , b = 1/3 implies u r ~ r 1 / 3 , which is Kolmogorov 
scaling. Here, we find that the moments of the entropy 
transfer rate (]Fe(k n )\ p ) also satisfy the same hierarchi- 
cal structure [22j and similarly r p is well approximated by 
Ci(l — 7 P ) — C2p. Similarly, C2 is the scaling exponent of 



lim r , 



{\F e (k n )\P +1 )/{\F (k n )\P). Follow- 



ing Ref. [20|, we estimate Fa (k n ) as S max u n k n , where 
Smax is the largest possible entropy. Since we observe 
BO-like scaling in the present case, it is more appropri- 

3/5 

ate to estimate u n ~ k n . As a result, we get C2 = 2/5. 
Then as (\Fg(k n )\) ~ (Fg(k n )) — x is independent of k n , 
we have t% = implying ci(l — 7) = C2. If we keep 
c\ = 1 as c = 1 for Xr, then we get 7 = 3/5. Putting 
these results together and using Eq. (fT4|) . we find 



C p -3p/5 = 1 - (3/5) p/5 - 2p/25 
?p-p/5 = 1 - (3/5) 2 ^ 5 - 4p/25 



(15) 
(16) 



Interesting, as shown in Fig. [2j Eqs. ([15]) and (fl~6|) indeed 
describe the measured intermittency corrections well. 

We have focussed on understanding the origin of 
anomalous scaling of an active scalar in homogenous tur- 
bulent convection using a dynamical shell model. We 
have extended Kraichnan's refined similarity idea to an 
active scalar in homogeneous turbulent convection and 
attributed the anomalous scaling to the variations in the 
entropy transfer rate. We have verified our hypothesis 
by showing explicitly that the conditional velocity and 
temperature structure functions at fixed values of the en- 
tropy transfer rate have simple scaling exponents of the 
BO values, and the intermittency corrections are given by 
the scaling exponents of the entropy transfer rate. Fur- 
thermore, by modifying earlier results obtained for the 
statistics of the local thermal dissipation rate in turbu- 
lent Rayleigh-Benard convection [20(, we have obtained 
the scaling exponents r p of the moments of the entropy 
transfer rate and thus Eqs. (TT5)) and <P~^|) for the intermit- 
tency corrections £ p — 3p/5 and £ p — p/5. These results 
are found to be in good agreement with the numerical 
values obtained in the simulations of the shell model. 

We should note that the scaling behavior of homoge- 
neous turbulent convection might not be the same as that 
in the central region of confined turbulent convection as 
coherent structures present in the latter case could af- 
fect the scaling properties [g]. Indeed direct numerical 
simulations [23j and analyses of experimental data 24 1 



indicated that the scaling behavior of the central region 
of confined turbulent convection is not well described by 
BO scaling plus intermittency corrections. On the other 
hand, there is evidence 25] of the validity of extending 
Kolmogorov's refined similarity idea in terms of the local 
thermal dissipation rate in the central region of confined 
turbulent convection. It would be interesting to further 
investigate this issue. 
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